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The observations of SNIa suggest that we live in the acceleration epoch when the densities of 
the cosmological constant term and matter are almost equal. This leads to the cosmic coincidence 
conundrum. As the explanation for this problem we propose the FRW model with dark matter and 
dark energy which interact each other exchanging energy. We show that the cubic correction to the 
Hubble law, measured by distant supernovae type la, probes this interaction. We demonstrate that 
influences between nonrelativistic matter and vacuum sectors are controlled by third and higher 
derivatives of the scale factor. As an example we consider flat decaying A(t) FRW cosmologies. We 
point out the possibility of measure of the energy transfer by the cubic and higher corrections to 
Hubble's law. The statistical analysis of SNIa data is used as an evidence of energy transfer. We 

\f} , find that there were the transfer from the dark energy sector to the dark matter one without any 

assumption about physics governing this process. We confront this hypothesis about the transfer 
with SNIa observations and find that the transfer the phantom and matter sector is admissible for 

£Nj , n m ,o = 0.27. We also demonstrate that it is possible to differentiate between the energy transfer 

model and the variable coefficient equation of state model. 
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PACS numbers: 98.80.Bp, 98.80.Cq, 11.25.- 



INTRODUCTION 

I i The modern cosmology, especially observational cosmology, reminds empirical science in a crisis phase Q. The 
\ recent observations of supernovae type la (SNIa) indicates that the Universe is accelerating at the present epoch 
0, El ■ We accept that the present evolution is well described by the general relativity theory with the Robertson- 
Walker type of space symmetry and the source of gravity is perfect fluid then the acceleration of Universe expansion 
can be explained only in the following way. The Universe is filled additionally to nonrelativistic matter with dark 
energy of unknown origin, which violates the strong energy condition. When the perfect fluid with the energy p and 
pressure p — the source of gravity — satisfies the strong energy condition then the explanation of SNIa observations 
requires the modification of the Einstein equations. If we postulate the Robertson- Walker symmetry then there are 
some propositions of modification of the Fricdmann first integral. Freese and Lewis considered flat cosmological 
models with the additional term Bp n , where B and n are constants The parameters of this model were confronted 
with the observation of distant type la supernova 0, @ . 

Both approaches can be tested statistically by searching model parameters which best fits to the SNIa data. But 
to distinguish among the models it is necessary to take additional observational constraints. When we assume the 
matter density value which is indicated by cosmic microwave background (CMB) and galactic counting observations 
then the most promising model is the dark energy model with the cosmological constant. The cosmological constant 
has very long history and is still the source of various problems and troubles 0, 13 ■ The main problem with today 
face of the cosmological constant is that its value is negligible in comparison to the Planck mass. In other words the 
A cold dark matter (ACDM) model looks like the effective theory which gives us the description of the phenomenon 
of acceleration without giving any understanding. 

Another problem is why the e nerg y densities of dark energy and of dust-like matter are of the same order of 
magnitude at the present epoch This problem is known as a "cosmic coincidence conundrum". The beginning 

of the vacuum energy and the dust-like matter is related to different epochs separated by very long interval of time. 
One of the possibility of the explanation of this coincidence is the intrinsic feedback between the energy density of 
dark matter and dark energy modelled by quintessence scalar fields ■ In the context of scalar quintessence fields it 
was proposed another interesting model with interaction It would be worthy to mention a model with coupling 
between dark energy and dark cold matter which reproduce power law solutions for energy density . This relation 
was constrained by the CMB observations Q| an d SNIa observations using statefinder diagnostic parameters (see 
|l5j and references therein). 

We differ in the presented approach that we do not assume any physical mechanism of the energy transfer, which 
is treated on the phenomenological level. We argue that luminosity distance versus redshift relation is the very first 
cosmological test that probes the interactions between dark matter and dark energy. This interaction is proposed and 
is checked whether the "cosmic coincidence conundrum" is solved. In particular we are interested in the direction of 
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energy transfer. For this aim we assume the Friedmann equation 

a 2 = ^a 2 -k (1) 
3 V ; 

where a is the scale factor, k is the curvature index, and a dot means the differentiation with respect to the cosmological 
time t. The second equation describing the evolution of the model is based on the adiabatic condition T£ = which 
for the Friedmann-Robertson- Walker (FRW) models with some perfect fluid assumes the form 

| (pa 3 )+p | (a 3 )=0 . (2) 

Because eq. |(5J| has the local character the standard interpretation is that p and p describe the effective energy density 
and pressure of multifluid which do not interact each other. Then eq. (J2J) describes separately the evolution of each 
component. 

Equation (J2J can be rewritten to the form 

P =-3H(p + p) (3) 

where the Hubble function H = a/ a. Let us note that the cosmological constant (for which p = —p) does not 
contribute to the conservation condition @ as long as it is treated as the noninteracting with the rest matter. 

We postulate that apart from dust matter there is dark energy, but both fluids interact now and energy can be 
transported from dark energy sector to the nonrelativistic matter sector. Therefore, eq. (jSJ cannot be separable for 
every component of multifluid. The special case of considered class of model are decaying vacuum cosmologies or A(t) 
models [lflll3|- The inspiration for constructing the noninteracting cosmologies is taken from Wojciulewitsch [13 (in 
context of dark energy see (lsjl. 

The main aim of the paper is to show that the observation of distant SNIa offer the possibility of testing the energy 
transport from the vacuum sector to the nonrelativistic matter sector which includes dark matter. We show that 
the measurements of third order term in the expansion of the luminosity distance relation with respect redshift z 
(jerk) allows to detect the energy transport. Higher order terms in the expansions (snap, crackle, etc.) control the 
velocity, acceleration of energy transport. Note that statefinder parameters also control third derivatives but they are 
inadequate if we want to detect the energy transfer directly from observations. 

We assume the different interpretation of eq. rather than its modification. 



TWO-SECTOR MODELS WITH TRANSFER OF ENERGY 



We construct the general class of the decaying dark energy models with the interaction starting from the Friedmann 
first integral which is independent of the form of pressure of fluid. We assume for simplicity some two-component 
fluid with effective pressure and energy 

PeS = PX + 0, Poff = p m + p x (4) 

where px — wxPx {wx = const) describes dark energy and p m is the energy of dust matter. If we put wx = — 1 the 
special case of the cosmological constant is recovered. 

The expression for the conservation condition can be rewritten to the form 

?l^« 3 ) + »l(«« 3(1+ " , )= < >- < 5 > 

The first term in eq. |5J) describes the net rate of absorption of energy per unit time in unit of comoving volume 
transfered out of the decaying vacuum fluid to the sector of nonrelativistic fluid. If we consider w = —4/3 the 
phantom fields are transported. Relation (JjJJ is usually interpreted without interaction between the sectors. Following 
Wojciulewitsch we postulate that the local energy conservation law JSJ can be written as 

^! =7(t) and ^^l(^ 3(1+ ^)=-^)- W 

The function j(t) is only a phenomcnological description of interaction between two sectors. Of course, the exact 
model of this interaction should be taken from the particle physics. If j(t) > the energy is transfered out of the 
vacuum, while if ~/(t) < the energy is transfered in the opposite direction. 
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Integration of eq. © gives 

Pm a 3 = Pm , a 3 + I 7 (t)a 3 di and p x a^ 1+w ^ = p x , a {1+Wx) - I 7 (t)a 3 < 1+ '°*><ft 



(7) 



where the index "0" means that the quantities are evaluated today. 

It would be useful for our further analysis to represent the Friedmann first integral Q in the form for a particle 
moving in the one-dimensional potential 
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Because of relation |(7J| the potential function is explicitly time dependent and now takes the following form 



V(a) 



A(t) , B{t) 



,1+3ujx 



(8) 



(9) 



where 



A(t) 



Bit) 



p x a 



3(l+w x ) 



In the concordance ACDM models both matter and the cosmological constant are treated separately without the 
interaction, so both functions A(t) and B(t) (densities p m ,o and px,o) are constant. The presence of the interaction 
manifests in the model by appearing the time dependence of the potential function JJjJ . 

Let us note that we postulate the time dependence of j(t) through the scale factor, i.e., j(t) = j(a(t)), and the 
potential function becomes only a function of a but the exact form of 7(a) is required. It is convenient to represent 
the dynamics of the model in terms of the Hubble function 



2 _ A(t) Bjt) 
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If we postulate that A(t) = A(a(t)), B(t) = B(a(t)) and put 1 + z = a -1 then relation can be used to fit the 
model parameters J7^o to the SNIa data, where i denotes all fluids considered. 

Differentiation of both sides of eq. JSJ for the potential in the form © gives the expression for acceleration 



a= 2 



A(a(t)) (l + 3w x )B(a(t)) 



2+3w x 



(11) 



where where we substitute derivatives A = ^j(t)a 3 and B = —^a 3 ( 1+Wx " l "f(t) = — a 3wx A from (|7|). Let us rewrite 
equation to the new form 



qH 2 



A(t) , (l + 3w x )B(a(t)) 



3(l+w x ) 



(12) 



where q = —-^1 is the deceleration parameter. 

To control higher derivatives of the scale factor we introduce the dimensionless parameter 



Qn = (-1 



,n+l 



1 d n a 
iH n dt n ' 



and then 



Q2 = q, 



3, 



(13) 



(14) 



are the deceleration, jerk, and snap, respectively [20||. In turn, to control the interaction we introduce the dimensionless 
transfer parameter 



3H 3 



(15) 



by analogy to the matter density parameter f2 m = p m /3H 2 
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To verify the model we estimate of the parameter v at the present epoch from the observation of SNIa data. From 
equations (|10|l and (|12(l we have 

„ k 3wy „ 2a , 

- Q k = = T~r-T — m + TlA 1 16 

a z H^ 1 + iwx 1 + owx 



and 



1 l + 3w Xn n „, 
g=-S2 m H il x - (17) 



Of course we have also the constraint condition fijt + i7 m + Qx = 1 • 

After the differentiation of both sides of eq. we obtain the basic equation relating the jerk to the transfer 
density parameter 

j - \w x v = 0„ + ^(1 + 3u;x)(2 + 3wx)fl x (18) 

and 

3 9 

j - ~wxv - 1 = -w x (w x + l)&x - fife. (19) 

Both for strings (wx = — 1/3) and topological defects (wx = —2/3) relation i|18fl does not depend on the density 
parameter of dark energy fix- This relation is obvious for all models. In the special case of the flat model fix = 1—Ct m 
and then we obtain 

3 9 

j- -w x v= -w x {w x + 1)Q X + 1. (20) 

Summing (|18|l and l|17|) for any we obtain 

3 3 3 

j - ^w x v + q= + jj(l + 3w x ){w x + 1)Q X (21) 

In the special case of the flat model formula l|21|l reduces to 

3 3 3 

j - WxV + q= --n m (4 + 3w x )w x + -(1 + 3w x )(l + w x ) (22) 

where wx can be always expressed in terms of fi m and q. Note that the relation 7(g) does not depend on priors on 
£! m for phantoms. 

Finally we obtain that the measurements of the jerk jo at the present epoch probes directly the effects of energy 
transfer as a consequence of relation l|21l) . While the cubic term in the relation cIl(z) is the first term in the Taylor 
expansion that depends explicitly on the "f(t), the higher terms in this expansion are related to the derivatives of "f{t). 
Let us define the parameter v n for the characterization of variability of j(t) as 

- - 1 (23) 



3#»+ 3 dt n 



It is obvious that vq = v, v\ = j/SH 4 , .... 

As an illustration that one can control the first derivative of ^(t) by the measurement of the snap Q4 = s, we prove 
the existence of some relation obtained by the differentiation of both sides of ill' 21) . To this aim we use the following 
formulas 

£ = + (24) 

' / --' n H[v + n m (2q-l)} (25) 



dt 
dn x 
dt 



-H[v + n x (l + 3wx -2q)}. (26) 
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Finally we obtain the relation 

3 3 3 

- s + j + 3jq + q + 2q 2 - -W X V\ = ^{3w x + Aw x ) + ^m(2q - 1) 

- ~(1 + 3w x )(l + w x )(l + 3w x - 2q)n x . (27) 

Let us briefly comment on the important case of wx = — 1 corresponding the decaying cosmological constant A(t) 
cosmologies. Of course, they constitute some special case of the considered models 



j + -v=l- fi fc 

■ 3 



(28) 
(29) 



The parameter j(t) = —dA/dt describes the first derivative of A and therefore the parameter v controls its variability. 
In turn the parameter v\ characterizes the convexity of the function A(t) 



3 9 3 

-s - j + 3jq + q + 2q 2 + -v x + 3v + -qv = -fl m (2q - 1) 



(30) 



TRANSFER PARAMETER FROM DISTANT SNIA 

Let us consider the luminosity distanced versus redshift relation di(^) expanded in the Taylor series with respect 
to redshift z. It can be done without knowledge about dynamical equation. For simplicity of presentation of the idea 
of measurement v, i>i, . . . we consider the flat universe what is justified by WMAP measurements. Then we obtain 

El 



4W = f 



1 + ~(1 - qa)z - i(l - q - 3gg + ] )z 2 



z 



24H 



[2 + 2q - I5q 2 



15q$ 



5j + lOgojo + s l 



(31) 

. .. Let us start with 
For this aim we mark 



We propose to detect the time variation of energy transfer using the parameters v, v\ , . 
estimation of v as a first approximation. We find the current constraints to the plane (go, J n) 
the shaded region of the 95% confidence level constraint from the recent SNIa measurements |22|. Because v = |(1— j) 
(or in the general case v = |(1 — j — £lk,o) the detection of the interaction is equivalent to the determination whether 
the jerk is different from 1 (or 1 — f2fe,o). If jo < 1 (or jo < 1 — f2fc,o) then the energy is transfered from the dark 
energy sector to the nonrelativistic matter sector. If jo > 1 (or jo > 1 — ^fc,o) the transport takes place in the opposite 
direction. Note that the negative curvature (fife,o > 0) makes the switch of transfer direction to happen for lower value 
of jo- Therefore, to find the direction of transfer we should know not only the value of jerk but also the curvature of 
space. 

In the general case for any wx we have 



3 3 

JO = -WX" + 2^m,0 



(32) 



where we use formula (II 7|) . 

We consider, for simplicity, the testing of the interaction for the flat model and the case wx = — 1 which corresponds 
the decaying cosmological constant. This allows to substitute v — > wxv. On Fig.Q] from relation H21(l with wx = — 1, 
the line jo = qo + §(^m.o ~ v ) is drawn when we assume that baryonic matter J7 mj o — v is equal 0.05. This relation 
allows us to estimate the interval on fl^o and jo on the la confidence level. We mark the line jo (go) & n d the vertical 
band to denote the interval with the la confidence level for qo e (—0.5725, —0.445) which gives fi m .o £ (0.285,0.37). 
In this interval of il m ,o the jerk jo is about 0.6. 

It is very interesting that present SNIa observations allowed us to measure the interaction without any special 
assumptions about physics of the transfer process. We thus determined the transfer energy parameter v and concluded 
that if we assume that the Universe is flat then the energy transfer takes place from the dark energy to dark matter. 
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n m0 = 0.285 
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•-• n m0 = 0.37 



100.0 
95.4 



O 



2 - 



- 



1 -' 




I 



68.3 



-1 




0.0 



-1.0 



-0.9 



-0.8 



-0.7 



-0.6 



-0.5 



-0.4 



-0.3 



-0.2 



-0.1 



0.0 



FIG. 1: The current constraint on the plane (go, jo)- The solid line is the relation jo(qo) fulfilled in the model. Additionally, 
for this relation the la confidence level interval for Q m ,o is drawn. 



In the previous section it was considered that Qb.o — ^m,o — v - Now we turn to estimation of the energy transfer 
parameter v using the SNIa data. It would be useful to consider two situation. First, the energy transfer is between 
decaying vacuum and matter sectors, and second, it is between the phantom (px = -i/3px) and matter sectors. We 
would like to answer on two questions. 

• What is the interval of Sl mj o which rules out the energy energy transfer [y = 0) on the confidence level 95%? 

• Is it possible to tell a scenario with energy transfer and another with variable u>x? 

To answer to the first question we test the hypothesis that v = 0. The transfer from decaying vacuum to matter 
sectors can be ruled out on the confidence level 95% for fl m .o 6 (0.23, 0.32) (Fig. [2J> , while the transfer from phantom 
to matter sectors is ruled out for f2 m ,o £ (0.30,0.37) on the same confidence level (Fig. [3J ■ Therefore the transfer 
between decaying vacuum and matter sector seems to be excluded because the extragalactic observations and CMB 
observations favor the values of f2 mi o in the obtained interval. On the other hand these other observations indicate 
that the transfer between phantom and matter sectors is possible. Hence if we have other arguments about phantom 
existence in the universe and if we accept that fl m .o ~ 0.27 as indicated by WMAP measurements then the energy 
transfer is necessary. 

Adopting the same analysis from the previous sections to the case of no energy and variable w(z) = wq + w\Z we 
obtain analogous formulas in which vwq is replaced by wittx.o- To answer the second question we analyze the Hubble 
diagram (Fig. It is shown that for very distant supernovae (z ~ 2) the model with variable w(z) predicts the 
brighter supernovae. 

The paper was supported by KBN grant no. 1 P03D 003 26. The author is very grateful dr A. Krawiec, dr W. 
Godlowski, and T. Stachowiak for comments and discussions during the seminar on observational cosmology. 
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